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The spectrum of Lyapunov exponents is powerful tool for the analysis of the complex system dynam- 
ics. In the general framework of nonlinear dynamical systems a number of the numerical technics 
have been developed to obtain the spectrum of Lyapunov exponents for the complex temporal be- 
havior of the systems with a few degree of freedom. Unfortunately, these methods can not apply 
directly to analysis of complex spatio-temporal dynamics in plasma devices which are characterized 
by the infinite phase space, since they are the spatially extended active media. In the present pa- 
per we propose the method for the calculation of the spectrum of the spatial Lyapunov exponents 
(SLEs) for the spatially extended beam-plasma systems. The calculation technique is applied to the 
analysis of chaotic spatio-temporal oscillations in three different beam-plasma model: (1) simple 
plasma Pierce diode, (2) coupled Pierce diodes, and (3) electron-wave system with backward elec- 
tromagnetic wave. We find an excellent agreement between the system dynamics and the behavior 
of the spectrum of the spatial Lyapunov exponents. Along with the proposed method, the possible 
problems of SLEs calculation are also discussed. It is shown that for the wide class of the spatially 
extended systems the set of quantities included in the system state for SLEs calculation can be 
reduced using the appropriate feature of the plasma systems. 
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Introduction 

Nonlinear dynamics and ehaos in plasma devices and 
microwave systems with interacted electron and electro- 
magnetic waves is often observed in bounded plasmas and 
electron and ions streams, where the nonlinearity and 
plasma instabilities are combined with dissipation caused 
by current flow to the boundary electrodes. Therefore, 
the development of the methods of analysis of spatially 
extended beam-plasma systems is a very important prob- 
lem due to such applications of chaos theory in plasma 
physics as the identification of chaoti spatial-temporal os- 
cillations in microwave and plasma devices (e.g. free elec- 
tron lasers, pasotrons, vircators, etc) [IHll, chaos control 
and suppression of turbulence [^-|8| , and chaotic synchro- 
nization in the plasma and microwave electronic systems 
ill. 

The Lyapunov exponents are powerful and useful tool 
for the complex system dynamics to be analyzed. They 
are used widely to characterize the systems with small 
number of deg rees of freedom in the different fields of 
science p^l - Uq . such as physics, physiology and medicine, 
molecular dynamics, astronomy, etc. 

At the same, the majority of the plasma systems are 
characterized by the infinite phase space, since they are 
the spatially extended active media. Unfortunately, the 
direct application of the Lyapunov exponent calculation 
technique developed for the system with the small num- 
ber of degrees of freedom to the objects of the beam- 
plasma and electron- wave systems being the spatially ex- 
tended systems is problematic. The main causes of the 
this problem are the following: (i) the phase space of 



the spatially distributed system is infinite; (ii) the num- 
ber of the Lyapunov exponents is also infinite; (iii) the 
system state (defined for each type of the system un- 
der study in its own way) must be used instead of the 
finite-dimensional vector; (iv) the orthogonalizaiton and 
normalization procedures based on Gram-Schmidt or- 
thonormalization must be modified to be applicable for 
the states. 

Nevertheless, the scientists try to use the Lyapunov 
exponents for the spatially extended systems due to the 
great efficiency of this tool. Mainly, all attempts are re- 
duced to the use of the modifications of the standard 
methods developed for the finite dimensional systems 
(e.g., the estimation of the highest Lyapunov exponent 
from time series or the representation of the spatially 
extended system as the finite dimensional one with the 
high dimension by means of the discretization). Unfor- 
tunately, these approaches do not take into account the 
distinguishing characteristics of the spatially extended 
systems and have serious limitations connected with the 
high dimension of the phase space of the discretizcd sys- 
tem 

Recently, the successful attempt to compute the high- 
est Lyapunov exponent taking into account the core fea- 
tures of the spatially extended systems has been made 
for plasma system models [l^ . In the present work we 
develop the technique of the spectrum of the spatial Lya- 
punov exponents (SLEs) calculation for the distributed 
beam-plasma systems and electron- wave devices. The 
structure of the article is the following. In SecUthe gen- 
eral idea of the proposed approach is given. The next 
Sec, mi contains the results concerning the computation 
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of the spectrum of spatial Lyapunov exponents for the 
spatially extended beam plasma system. As the model 
system under study the Pierce diode has been selected, 
with the spectrum of Lyapunov exponents being calcu- 
lated both for autonomous oscillator and for two coupled 
plasma diodes. To illustrate the generality and efficiency 
of the proposed method, in Sec. Illll the SLEs have been 
also calculated for the transverse field backward wave os- 
cillator being the classical models of the electrons-wave 
systems. The problem concerning the selection of the 
system state for the SLEs computation is discussed in 
Sec. IIVI The final remarks are given in Conclusion. 

I. THEORY OF THE SPATIAL LYAPUNOV 
EXPONENTS 

The main idea of the computation of the spectrum of 
the spatial Lyapunov exponents is the following. Let 

L(U(x,i))=0 (1) 

be an operator determining the evolution of the spatially 
extended system under study, \J{x,t) is the system state 
which (in the general case) depends on the spatial coor- 
dinate X G [0,L] and time t. 

The state U(a;, t) plays the role of the reference point 
to study the evolution of the small perturbations. To 
compute N largest spatial Lyapunov exponents we con- 
sider the set of perturbations Vi{x,t), i = I, . . . , N which 
are orthogonal 

(v.v,).{J; (2) 

where (U, V) is the scalar product. If the system state 
is the A'ti-dimensional vector U = (ui, . . . , ujv^)^ or real 
function U = u{x,t) defined in the definition of the 
scalar product is 

(U,V)=^^.,(i)^^.W (3) 

i=l 

or 

(u.v). (4) 

L 

respectively. At the same time, the state of the spatially 
extended beam-plasma system should be usually defined 
in a more complicated way, and, as a consequence, the 
definition of the scalar product is also more complex. 
E.g., if the system state is the iVrf-dimensional vector- 
function XJ{x,t) = {ui{x,t), . . . ,UNa{x,t))'^ the scalar 
product should be defined as 




Additionally to the orthogonality requirement ([2]) all 
perturbations Vi(a;,t) must also obey the normalization 
condition 

l|V,(a;,0)|| = l, (6) 

where ||V|| = v^(V,V). 

The set of the perturbations 'Vi{x,t) fulfilling require- 
ments ([2]) and (|n|) may be obtained with the help of the 
Gramm-Schmidt orthonormalization method 

' ^ V.(x,0) ^ V,{x,0) 

0V,(a:,O),V,(a:,O)) 

< Yi{x,0) = ^i{x) 

i 

V,+i{x,0) = </p,+i(.t) - ^(Vfc,v?,+i)Vfe(x,0), 
fe=i 

(7) 

i = 1, 2, . . . , iV — 1, where ipi{x), ip2{x), . . . , ^n{x) is the 
set of the linearly independent arbitrary states defined for 
the system under study. 

To compute the spectrum of the spatial Lyapunov ex- 
ponents of the spatially extended system the behavior 
both of the state U(x, t) of the system under study and 
all perturbations 'Vi{x, t) {i = 1,2, . . . , N — 1) has to be 
considered. The evolution of the system state U(a;, t) is 
defined by Eq. ([T]) , whereas the development of the per- 
turbations Vi{x,t) of the reference state is described by 
the linearization 

dL{V{x,t),V,{x,t)) ^0 (8) 

of the evolution operator ([T]) applied in the vicinity of 
the state U{x, t). 

After the time interval with the length T the obtained 
set of the perturbations Vi(a;,T) must be orthogonal- 
ized and normalized with the help of the Gramm-Shmidt 
method again, with the set of the functions ifi{x) being 
the considered perturbations Vi(a;,T). In other words, 

if,{x)=Y,{x,T). (9) 

The described algorithm for spatially extended system 
must be repeated many times. After M iterations the 
Lyapunov sums are calculated 

M 

5, ==^ln||V,(x-,jT)||, (10) 
i=i 

where the perturbations before renormalization but after 
orthogonalizaiton are used. The spatial Lyapunov expo- 
nents are estimated as 

A, = — . (11) 

II. SPECTRUM OF SPATIAL LYAPUNOV 
EXPONENTS FOR PIERCE DIODE 

To illustrate the efficiency of the proposed method for 
the computation of spatial Lyapunov exponents we con- 
sider the behavior of autonomous and coupled Pierce 
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FIG. 1: (Color online) Schematic diagram of Pierce diode 



diodes [17[ being the classical models of beam-plasma 
systems, demonstrating the complex spatio-temporal os- 
cillations including the chaotic ones [l7l - [l9j . 



Pierce diode 



A Pierce diode |17H20l | is schematically illustrated in 
Fig-Hl and consists of two infinite parallel plains pierced 
by a mono-energetic electron beam. Grids arc grounded, 
with the distance between them being L. The charge 
density po and electron velocity dq are constant at the 
system input. The region between two plains is uniformly 
filled by neutralizing stationary ions, with the density 
IPionI being equal to the non-perturbed electron beam 
density |po|- 

The dimensionless Pierce parameter a ~ ujpL/vQ de- 
termines the dynamics of the system (here ujp is the elec- 
tron beam plasma frequency, vq is the non-perturbed 
electron velocity, L is the distance between the diode 
plains). With a > n the so-called Pierce instability is 
developed in the system and a virtual cathode is formed 
in the electron beam [l^. A detailed analysis of the 
Pierce instability is provided in Ref. [21 1. 

Here we focus our attention on a narrow range of Pierce 
parameter near a ~ Stt where the increase of the instabil- 
ity is suppressed by the nonlinearity, and a reflectionless 
regime takes place in the electron beam [l^, [2^ . In this 
case the s yste m behavior may be described by the fluid 
equations |18l - [2g | , that give rise to various types of beam- 
plasma chaotic oscillations 0, 0, [l^, [22| ■ 

The dynamics of the electron flow in the Pierce diode 
(in the fluid electronic approximation) is described by the 
self-congruent system of dimensionless Poisson, continu- 
ity and momentum equations 



ay 



dp dp dv ^ 

dt ^ dx ^ dx 



dv dv 
■ ■ V- 



d(p 
dx ' 



(12) 
(13) 
(14) 



dt dx 
with the boundary conditions 

viO,t) = l, p(0,i) = l, (^(0,t) = (p(l,t) =0, (15) 



where (p{x, t) is the dimensionless potential of the electric 
field, p{x, t) and w(.t, t) are the dimensionless density and 
velocity of the electron beam (0 < x < 1), respectively. 
The dimensional variables ((p', p' , v' , x', t') are connected 
with the dimensionless ones as <p' = (y'^/ri)ip, p' = pop, 
v' = wow, x' = Lx, t' = {L/vo)t, where 77 is the specific 
electron charge, vo and po are the non-perturbed velocity 
and density of the electron beam. 

The system (ll2p - (fT4l) with the boundary condi- 
tions (Uni) describes the dynamics of electron fiow in the 
autonomous Pierce diode. It has been solved numeri- 
cally using finite difference approximation. Numerical 
solutions for motion ([T4| and continuity (fT3|) equations 
have been found by means of the explicit scheme with 
the differences against flow, and Poisson equation (fT2|) 
has been integrated using the error vector propagation 
method [2^. The time and space integration steps have 
been selected as Ax = 0.005 and At = 0.003, respec- 
tively. 

As a vector state of Pierce diode we have used the 
following one 



U= (p(a;,t),w(x,t))^ 



(16) 



In this case, according to ([5]), the scalar product is defined 
as 

1 

(Ui,U2) = J[pi{x,t)p2{x,t)+vi{x,t)v2{x,t)]dx, (17) 



where the integration in (jl7p should be performed over a 
whole length of the considered system. Note, that the po- 
tential ip{x, t) is not included in the vector state U since it 
could be explicitly defined through the density of the elec- 
tron beam p{x, t) by Poisson equation (jl2p and boundary 
conditions In other words, to explicitly define the 

Pierce diode system state it is sufficient to know the dis- 
tributions of the density p{x^ t) of the electron beam and 
velocity u(x, t). It allows us, using the evolution operator 
i(U), i.e. the system ([T^ - p^ with boundary conditions 
([T5|). to calculate the analogues states in the other mo- 
ments of time. 

As the perturbations we have used the vectors 

(cf(x,t),en^,i))'^, 



(18) 



characterizing the small deviations from the reference 
state U, with the small perturbation £^f{x^t) of the po- 
tential <p(a;, t) being explicitly defined by the distribution 
of the density of the electron beam p{x,t) and its small 
perturbation ^f(a;,t). 

The behavior of Pierce diode with i-th perturbation of 
the state U is described by the following system 



with the boundary conditions 







p(0,t)+C''(0,t) = l, 

^(0,*)+^" (0,0 = 1, 
(p(0,t)+Cf(0,t) = 0, 

^(1,0+^(1,0 = 0. 



(19) 



(20) 



4 



Having linearized the system (jl9[) we can rewrite the 
equations ([T^. (PU)) in the fohowing way 



aL(u,v,) = 0, 



(21) 



where 9L(U, Vj) is the evolution operator L(-) linearized 
near the reference state \J{x,t). It describes the char- 
acter of the development of the perturbations with 
time. For the system ((T2|) - (fT4)) the evolution operator 
dL{lJ,Yi) can be written in the following way 



3cV 



- CP 



dt ^' dx 



dx 



dx 



dx ' 



(22) 



dt dx * dx dx 

In this case the boundary conditions (j20p for the small 
deviations of the synchronous state are given by 



C(o,<) = o, 

er(o,o = o, 



(23) 



Obtained operator dL{\J, V) is linear to relatively small 
deviations V. 

Now, taking into account the reference state U, the 
perturbations Vi, the evolution operator L(\J) defined 
for system p^ ~ ([T4)) with boundary conditions ([T5|). 
and the linearization operator 9L(U,V) given by rela- 
tions we can calculate the spectrum of spatial 
Lyapunov exponents for the autonomous model of Pierce 
diode. 

The behavior of autonomous model of Pierce diode is 
known to be defined by the only one dimensionless con- 
trol parameter a called Pierce parameter. As we have 
mentioned above, Pierce diode is capable to demonstrate 
the different types of the complex spatio-temporal oscil- 
lations, with the transition from the periodic to chaotic 
ones being performed via t he p eriod doubling bifurca- 
tion cascade [H, [l^ [H, [l^, H^] when the a-parameter 
value decreases in the range a G (2.8627r, 2.887r). Fur- 
ther decrease of the parameter a results in the consid- 
erable complication of the chaotic regime being realized 
in system ([T2 |) - (fT4l) attended by the disappearance of the 
well-defined time scale of the spatio-temporal oscillations 
of the electron beam with the further transition to the 
wide-band chaotic spatio-temporal oscillations. 

To characterize the type of dynamical regime being 
realized in the system under study we compute the spec- 
trum of spatial Lyapunov exponents. In Fig. [51 a the 
dependence of the four highest spatial Lyapunov expo- 
nents on the parameter a is shown. To calculate the val- 
ues of the spectrum of the spatial Lyapunov exponents 
M = 2 X 10^ iterations of the renormalization and orthog- 
onalizaiton have been done for each value of the control 
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FIG. 2: (Color online) (a) The dependence of the four high- 
est Lyapunov exponents and (6) bifurcation diagram of os- 
cillations in the autonomous model of Pierce diode with a- 
parameter value increasing 



parameter a, with the time interval between iterations 
being equal to T = 0.6. In parallel, the bifurcation dia- 
gram of the density of the electron beam oscillations at 
the point of the interaction space x = 0.2 is also indicated 
in Fig. [2J&. It is easily seen from Fig. [5] that the spec- 
trum of the highest spatial Lyapunov exponents allows to 
detect clearly the type of the dynamical regime realized 
in the system under study. For all values of the control 
parameter a one zero Lyapunov exponent corresponding 
to the system perturbations along the time axis is pre- 
sented in the spectrum. In the region of the periodic os- 
cillations (e.g., a < 2.862) such zero Lyapunov exponent 
is the highest one, i.e., the chaotic dynamics is absent. 
All other spatial Lyapunov exponents are negative. At 
the same time, it is clearly seen that in the bifurcation 
points corresponding to the period doubling cascade (see, 
e.g., a ~ 2.87) the second Lyapunov exponent A2 is also 
equal to zero, that indicates the qualitative changes in 
the system behavior with the control parameter varying 
as well as for the low-dimensional dynamical systems. 

With the further decrease of the parameter a (see 
Fig. [2]) the highest spatial Lyapunov exponent Ai be- 
comes positive, i.e. the transition to chaotic oscillations 
takes place. The second Lyapunov exponent remains to 
be the zero one. All others Lyapunov exponents are neg- 
ative. Note, that in the periodicity windows, where the 
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periodic oscillations take place, the positive Lyapunov 
exponent becomes also negative and the zero Lyapunov 
exponent gets to be the highest one. 

It is important to note that in the case of the chaotic 
oscillations the system under study has the only one pos- 
itive Lyapunov exponent. At the same time, the spa- 
tially extended media may theoretically have several pos- 
itive Lyapunov exponents. The presence of the only one 
positive Lyapunov exponent is the evidence of the sim- 
plicity of chaotic oscillations being realized in the spa- 
tially extended system. One can assume that just due 
to that fact it has been possible to construct the finite- 
dimensional model of the dynamics of the electron beam 
in the Pierce diode being the system of ordinary differ- 
ential equations [lO] by means of Galerkin method (23l |. 

So, the spatial Lyapunov exponents could be easily ap- 
plied for the detection and classification of the dynamical 
regimes being realized in the beam-plasma systems with 
the overcritical current. However, the field of the possi- 
ble applications of the spectrum of the spatial Lyapunov 
exponents is not restricted to the consideration of the au- 
tonomous system. In particular, Lyapunov exponents are 
an effective tool for the analysis of the different types of 
chaotic synchronization being realized both in the finite- 
and infinite-dimensional systems [26l - [28| . Therefore, we 
consider the spectrum of spatial Lyapunov exponents for 
two unidirectionally coupled Pierce diodes in the next 
subsection. 



B. Two coupled Pierce diodes 

The system under study in this case is given by the 
partial differential equations 



dt 



dx 



,dv^ 



dip^ 



dx 



(24) 
(25) 
(26) 



dt dx 
with the boundary conditions 

t;i-2(0,<) = l, p''\0,t)^l, /•2(0,t)=0, (27) 

where the top indexes "1" and "2" correspond to the first 
and second coupled beam-plasma systems, respectively. 

The unidirectional coupling between considered Pierce 
diodes is realized by the modification of the boundary 
conditions on the right boundary of the second (response) 
system, exactly in the same way as it has been done 
in In other words, the boundary conditions for the 
system ((24)) - ([26)) describing the dynamics of unidirec- 
tionally coupled Pierce diodes should be written as 



ip\l,t) = 

(^2(1, t) = £{p^X^l,t)~ p\x^l,t)), 



i.e., the first (drive "1") system is in the regime of au- 
tonomous oscillations influencing on the second (response 
"2") one. 

To analyze the behavior of two coupled Pierce diodes 
we compute spatial Lyapunov exponents for the system 
(|24| - (|26)) . As a system state we use now the vector 



(29) 



U = {p\x, t), v\x, t), p^{x, t), v^{x, t)f 



According to Eq. ([5]) the scalar product of the system 
states in this case is defined as 

1 

(Ui,U2) = / {p\{x,t)p\{x,t) + v\{x,t)vl{x,t)) dx+ 



1 

/ {pl{x,t)pl{x,t)+vf{x,t)v^{x,t)) dx. 



(30) 

In Eq. (|30)) the integration is performed on the whole 
length of the considered system as well as in the case of 
Eq. (H?)). Equations ((MJ-din]) with the boundary con- 
ditions (|27)) - (|28| for two unidirectionally coupled Pierce 
diodes play the role of the evolution operator ([T|) . 
By analogy with the relation (fT8)) we use vectors 

V, = i^l''{x,t),^mx,t),^^''{x,t),^f^x,t)f, (31) 

as the perturbations for two unidirectionally coupled 
Pierce diodes. They characterize the small deviations 
from the reference state U, with the small perturbations 
^l'''{x,t) and ^'f'^{x,t) of the potentials ip^''^{x,t) as well 
as in the case of autonomous Pierce diode being explicitly 
defined by the distributions of the density p^''^{x, t) of the 
electron beam and its small perturbations {x^t). 

For the system under study ()24p -((28 )) the operator 
(9L(U, V;) defining the evolution of the small perturba- 
tions V; of the reference state U is written in the form 



dx"^ 




1 




dt 


fl,2pdv''' 

dx 

A,2v V'^ 

dx 


qA,2p 

OX 
1,2 ■ 

^ dx ' 




dt 


OX 


A,2v 9v^'^ 
dx 


dx 



(32) 



with the boundary conditions for the small deviations 
from the reference state being the following 



^'■''(0,0 = 0, 
e!'"'(o,t). = o, 

^!''^(0,i) = 0, 
^""(1,0-0, 

e^(l,t) = e(e(l,t)-e(l,i))- 



(33) 



(28) 



So, the computations of the spectrum of the spa- 
tial Lyapunov exponents for two unidirectionally coupled 
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FIG. 3: (Color online) The dependence of the five highest 
Lyapunov exponents on the coupling parameter e for two uni- 
directionally coupled Pierce diodes 



Pierce diodes is performed in the same way as for the au- 
tonomous system with the only one difference consisting 
in the fact that the system state U, small perturbations 
Vi, the scalar product (Ui, U2) and the evolution opera- 
tors L{XJ) and 9L(U, V^) are defined for two interacting 
systems. 

The dependence of the five highest spatial Lyapunov 
exponents on the coupling parameter strength e for two 
unidirectionally coupled Pierce diodes with the 

boundary conditions (|27l) ~ ([28|) is shown in Fig. [S] The 
control parameter values = 2.858, = 2.860 have 
been chosen in such a way that the autonomous dynamics 
of each Pierce diode is being chaotic one. It is clearly seen 
that in the absence of the coupling e = two highest Lya- 
punov exponents are positive (that is the evidence of the 
presence of the chaotic dynamics in each Pierce diode), 
two Lyapunov exponents are equal to zero (that corre- 
sponds to the perturbations of the system states along 
the time axis), the fifth and others Lyapunov exponents 
are negative. 

With the coupling parameter value increase the quan- 
titative values of the Lyapunov exponents start chang- 
ing. Due to the unidirectional type of the coupling be- 
tween Pierce diodes the spectrum of Lyapunov exponents 
Ai > A2 > . . . can be divided into two parts. One of 
them corresponds to the Lyapunov exponents calculated 
for the drive system Aj > > . . . and due to the inde- 
pendency of the behavior of the drive system on the re- 
sponse one they do not depend on the coupling parameter 
strength. The second part is the conditional Lyapunov 
exponents, i.e., Lyapunov exponents computed for the re- 
sponse system. They depend on the coupling parameter 
strength e and, therefore, are called conditional. 

Such a model of Pierce diode is known to demonstrate 
the different regimes of chaotic synchronization, includ- 
ing complete, generalized and time scale synchroniza- 
tion [1^ . The occurrence of any type of the synchronous 
regime should be accompanied by the transformations 
in the spectrum of Lyapunov exponents. In particular. 



for the systems with a small number of degrees of free- 
dom the transition of the zero conditional (corresponding 
to the response system) Lyapunov exponent in the field 
of the negative values precedes the phase or time scale 
synchronization regimes in the case of the small values 
of the parameter detuning [28l - [30j . When the highest 
conditional Lyapunov exponent passes through zero, the 
generalized synchronization regime takes place (l2l . Isij . 

The very same phenomena is observed for two cou- 
pled spatially extended beam plasma systems. When 
the coupling parameter value increases the response sys- 
tem becomes synchronized by the drive one, with the 
qualitative changes in the spectrum of the conditional 
Lyapunov exponents taking place. In Fig. [3] two crit- 
ical points £0 ~ 0.01 and e+ « 0.09 are clearly seen. 
As in the case of the low-dimensional system, the first of 
them corresponds to the transition of the zero conditional 
Lyapunov exponent in the field of the negative values, 
whereas the second one relates to the zero-pass of the 
positive conditional Lyapunov exponent. The negative- 
ness of the largest conditional spatial Lyapunov exponent 
is the evidence of the generalized synchronization regime 
that has been confirmed with the help of the auxiliary 
system method [H, [11] . 



III. SPECTRUM OF SPATIAL LYAPUNOV 
EXPONENTS FOR THE TRANSVERSE FIELD 
BACKWARD WAVE OSCILLATOR 

To illustrate the generality and efficiency of the pro- 
posed method, as well as their applicability for the 
widespread class of spatially-extended systems, we con- 
sider the behavior of the transverse field backward wave 
oscillator being the classical models of the electrons- wave 
systems, demonstrating the complex oscillations includ- 
ing the chaotic regimes as well as hyperchaos. 

The dynamics of the backward wave oscillator is de- 
scribed by the system of the dimensionless equations 



dF dF 
dt dx 



'AF 



dx 



+ J\I 



-AF 



(34) 



(35) 



where F(x, t) is the dimensionless complex- valued am- 
plitude of the electric field and I{x, t) is the complex- 
valued current (0 < a; < 1). The control parameter A 
can be considered as the dimensionless length of the sys- 
tem. The equations (|34 |) - (|35l) have been integrated with 
the boundary conditions: 



F(l,t) = 0, /(0,t) = 0. 



(36) 



To calculate the spectrum of Lyapunov exponents we 
define the main state of this system as 



U = F{x,t). 



(37) 
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It is important to note, that the distribution of the 
complex-valued current is excluded from the main sys- 
tem state, since it is explicitly defined by the electrical 
field F{x,t). 

The scalar product of complex- valued functions in this 
case is defined as follows 



(Ul,U2) 



Fi{x,t)F^{x,t) dx, 



(38) 



where the integration, as for the case of Pierce diode, 
should be performed over the whole length of the consid- 
ered system. 

The main system state (|37|) evaluates according to the 
operator ([T]), presented in the form of differential equa- 
tions (p4|) -([35 l) . These equations (p4|) - ([35|) have been in- 
tegrated using the Lax-Wendroff method and the third- 
order Runge-Kutta method. The time and space in- 
tegration steps have been selected as Ax = 0.004 and 
At = 0.002, respectively. 

The state (|57|) plays the role of the reference state. 
To calculate N largest spatial Lyapunov exponents we 
consider the set of perturbations 



V, = (ef(x,t)), 



(39) 



Again, the current has been eliminated from the con- 
sideration, since the small perturbations (a;, t) of the 
complex- valued current I{x,t) are explicitly defined by 
the distribution of the high-frequency electric field F{x, t) 
and its small perturbations {x,t). 

The evolution of the set of small perturbations is de- 
scribed by the linearized operator ([5]) , which for the con- 
sidered system can be written in the form: 



dtj dx ' 

^+2j\i\ur+ji'e.=-Aif 

with the boundary conditions 

F(l,t) + ef(l,t)-0, /(0,t) + e/(0,t) = 0. 



(40) 



(41) 



The Lyapunov sums and values of SLEs have been esti- 
mated with the help of expressions plI)) - ([TT]) . 

In Fig. 131 a the dependence of the four highest spatial 
Lyapunov exponents on the parameter A is shown. To 
give the information of the behavior of the lower SLEs the 
dependence of the 10th LE is also given in Fig. 01 a. The 
bifurcation diagram of the complex-valued amplitude of 
the electric field is presented in Fig. |4l6. The largest 
Lyapunov exponent Ai is equal to zero when the control 
parameter has the value being less then 1.83, that corre- 
sponds to the stationary state, when the antiphase oscil- 
lations of the real and imaginary parts of the complex- 
valued amplitude of electric field takes place, while its 
module remains constant. All other SLEs A2, A3, . . . are 
negative. 



A 
2 



-2 

-4 



A,' 




■ JfAio 





1.5 



2.0 2.5 3.0 3.5 4.0 4.5 




1.5 2.0 2.5 3.0 3.5 4.0 4.5 A 



FIG. 4: (Color online) (a) The dependence of the four highest 
and the 10th Lyapunov exponents (the second LE is shown 
by points (•)) and (b) the bifurcation diagram of oscillations 
in the transverse field Backward wave oscillator with the A- 
parameter value increase 



For the bifurcation point Ac « 1.83 the stationary 
state becomes unstable and, as a consequence, the pe- 
riodic oscillations of the module of the complex-valued 
electric field \F\ are observed in the system under study. 
In the region of the periodic oscillations A > Ac the sec- 
ond Lyapunov exponent A2 is also equal to zero. At the 
same time, all other Lyapunov exponents are negative, as 
before. Note also, at the bifurcation point Ac the third 
Lyapunov exponent A3 is equal to zero, that indicates the 
qualitative changes taking place in the spatially-extended 
system. 

With the further increase of the A-parameter the tran- 
sition from periodic oscillations to chaotic dynamics is 
observed. This transformation is accompanied by the 
modification of the spectrum of SLEs, when the largest 
Lyapunov exponent becomes positive {Ach ~ 3.9). As 
for the case of Pierce diode, it can be used as a criterium 
of the transition to the chaotic regime, and, as well as for 
Pierce diode, this transformation of the SLEs spectrum 
agrees with the bifurcation diagram (Fig. |31a). 

As it has been described above, the chaotic regimes in 
the spatially-extendent systems, in general, can be char- 
acterized by several positive Lyapunov exponents. It is 
clearly seen from Fig.|3l6, the number of the positive Lya- 
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punov exponents increases with the grow of the value of 
A-parameter, with the chaotic regime (see bifurcation di- 
agram, Fig. IHa) becoming more comphcated. Note also, 
as soon as one of the Lyapunov exponents being zero be- 
comes positive, the largest negative Lyapunov exponent 
takes the value zero. 



IV. THE SELECTION OF THE SYSTEM STATE: 
BRIEF DISCUSSION 

The proposed approach for the computation of the 
spectrum of the spatial Lyapunov exponents is shown 
to be the effective and powerful tool to study and ana- 
lyze the complex behavior of the beam plasma systems 
as well as the other spatially extended media. Therefore, 
it is not surprised that this tool attracts great attention 
of scientists [H, [sj]. Nevertheless, one have to be very 
careful computing the spectrum of the spatial Lyapunov 
exponents, since there may be certain particularities con- 
nected with the SLEs calculation which may result to the 
inaccurate or even incorrect decisions made and, there- 
fore, they must be taken into account by scientists. In 
this section one of the possible problems will be discussed 
for the scientists to avoid the potential errors. 

When the spectrum of spatial Lyapunov exponents is 
calculated the state U of the system should be chosen 
properly. In the studies given above certain variables 
(such as the potential (p{x,t) in the case of Pierce diode 
or the current I{x,t) in the case of the transverse field 
backward wave oscillator) have been excluded from the 
system state having based on the reason that these vari- 
ables may be explicitly defined through the other ones 
which are included in the system state. To prove the va- 
lidity of this approach we have also calculated the spec- 
trum of the spatial Lyapunov exponents when the whole 
set of the variables is supposed to be the system state. 
For the Pierce diode the system state in this case is 

IJ = {fix, t),p{x,t),v{x,t)f, (42) 

whereas for the transverse field backward wave oscillator 
it is defined as 

lJ^iF{x,t),I{x,t)f. (43) 

Fig. [5] illustrates the three highest Lyapunov expo- 
nents of Pierce diode obtained in both cases (when U = 
{(p{x,t),p{x,t),v{x,t))'^ and U = {p{x,t),v{x,t))'^). 
One can see that the LEs correspond exactly to the val- 
ues, calculated in Sec.|lll Analogously, in Fig. [6] the first, 
the third and the fifth Lyapunov exponents correspond- 
ing to the periodical oscillations in the backward wave 
oscillator are shown. It is clearly seen, that the val- 
ues of Lyapunov exponents obtained in the both cases 
(U = {F{x,t),I{x,t))'^ and U = F{x,t)) are identical 
for all values of the dimensionless system length A. 

So, certain variables can be excluded from the sys- 
tem state if they may be explicitly defined through the 



2.855 2.865 2.875 2.885 a 



FIG. 5; (Color online) The three highest Lya- 
punov exponents of Pierce diode calculated both 
for U = {ip{x,t), p{x,t),v{x,t))'^ (points) and for 
U — (p(x,t),v{x,t))'^ (lines) 

other system variables in the every moment of time or, 
in other words, if the time derivatives of these variables 
are not presented in the evolution operator. Neverthe- 
less, in general, the choice of the system state U requires 
the thoroughness and accuracy. Thus, the majority of 
the beam plasma and microwave systems are simulated 
with the help of the particle methods [s^. In this case 
the state U(x, t) of the system under study is determined 
by the electromagnetic field as well as locations and ve- 
locities of all particles used for simulation. As far as the 
electromagnetic field is concerned there is no problem 
to apply the rcnormalization and orthogonalizaiton pro- 
cedures in the very same way as it has been described 
in Sec. HI since this field is the function of the coordi- 
nate, i.e., E{x,t) € L^, B{x,t) € L^. On the contrary, 
the set of particles emulating the electron beam can not 
be renormalized and orthogonalized with the help of the 
technique given above, since each particle is characterized 
by its own coordinate and velocity. As a consequence, the 
scientists dealing with such a type of systems use only the 
electromagnetic field as the system state U(a;, t) to avoid 
the uncertainty connected with particle dynamics eval- 
uation (see, e.g. H^). Although in some instances the 
obtained in such a way results seem to be plausible, we 
show below that this approach gives the very crude and 
inaccurate solution. 

To illustrate this point we consider the emulation of the 
spatial Lyapunov exponent calculation for the spatially 
extended system simulated by the particle method by 
the example of Pierce diode. Since in this case the elec- 
tron beam ( "particles" ) can not be used when the rcnor- 
malization and orthogonalizaiton procedures are imple- 
mented, wc consider the truncated system state 

V'{x,t) ^v{x,t) (44) 

in parallel with The dynamics of Pierce diode state 
and its perturbations are evaluated according to the op- 
erator P^ - p4[) . whereas the Gramm-Schmidt procedure 
is realized only for the truncated state pij) that excludes 
the beam of electrons from the rcnormalization and or- 
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FIG. 6: (Color online) The first, the third and the fifth 
liighest Lyapunov exponents of the backward wave oscilla- 
tor calculated both for U — {F {x , t) , I {x , t))'^ (points) and 
U = F{x, t) (lines) 
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FIG. 7: (Color online) The four highest Lyapunov exponents 
calculated for the case when the beam of electrons is ex- 
cluded from the renormalization and orthogonalizaiton pro- 
cedure (compare with Fig. [2] a) 

thogonalizaiton, exactly in the same way, as it has been 
done, e.g., in [33| for the system simulated by the particle 
method. 

The four largest Lyapunov exponents (solid lines) cal- 
culated for Pierce diode by means of the described tech- 
nique (it models the approach used for the system sim- 
ulated by particle method) arc shown in Fig. [7] to be 
compared with the Lyapunov exponents calculated in 
Sec. Ill Al (see Fig. [5]). It is easy to see that the true 
dependence of the spatial Lyapunov exponents (Fig. [2]) 
is seriously distorted (Fig. [T]) when the electron beam is 
excluded from the renormalization and orthogonalizaiton 
procedures. E.g., even zero Lyapunov exponent being 
known to have always to exist in the autonomous system 
(since it corresponds to the perturbation along with the 
time axis) is not calculated properly in comparison with 
the case, when all sufficient quantities (in the considered 



case these are p{x,t) and v{x,t)) are taken into account. 

Evidently, for some values of the control parameter a 
the obtained set of spatial Lyapunov exponents seems to 
be reasonable. E.g., for a = 2.856 the largest value of the 
calculated spatial Lyapunov exponent spectrum is posi- 
tive that corresponds to the chaotic dynamics observed 
in the system, the second Lyapunov exponent is close 
to zero that also agrees well with the theory, the other 
Lyapunov exponents are negative. So, if the careful ver- 
ification of the obtained results has not been done these 
results may be assumed to be correct, although the val- 
ues of the calculated spatial Lyapunov exponents differ 
greatly from the true ones. Moreover, for the other values 
of the control parameters the calculated data may be ob- 
viously contradictory (see, e.g., a = 2.87, where despite 
of the periodic regime the positive value of the spatial 
Lyapunov exponent is obtained, whereas the zero Lya- 
punov exponent is absent). The inconsistency of the ob- 
tained values of the spatial Lyapunov exponents with the 
system dynamics may be revealed with the help of the 
examination of the dependence of SLEs on the system 
control parameter value, but, typically, this analysis is 
not given. 

Therefore, we can come to the conclusion that the cal- 
culation of the spatial Lyapunov exponents for the mod- 
els simulated by the particle methods must be used very 
carefully. Unfortunately, at present there is no methods 
allowing the calculations of SLEs for such kind of the 
dynamical systems and this problem should be solved in 
the future. 



Conclusion 

In the present work we have proposed the method for 
the calculation of the spectrum of the spatial Lyapunov 
exponents for the spatially extended beam plasma sys- 
tems. This method is shown to be effective and powerful 
tool for the analysis of the complex behavior of the beam 
plasma systems as well as the other spatially extended 
systems. The possible applications of the spectrum of 
the spatial Lyapunov exponents are not restricted to the 
consideration of the autonomous system, e.g., SLEs are 
the effective tool for the analysis of the different types 
of chaotic synchronization. At the same time, the cal- 
culation of the spatial Lyapunov exponents for the mod- 
els simulated by the particle methods must be used very 
carefully, whereas the problem of calculations of SLEs for 
these systems should be solved in the future. 

Though the calculation of the spectrum of the spatial 
Lyapunov exponents has been illustrated with the help of 
several examples, we expect that this technique may be 
used in many other relevant circumstances, as e.g. laser 
systems [11], semiconductor superlattices [s^l, etc. 
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